Dense Colloidal Suspensions under Time-Dependent Shear 



J. M. Bradcr\ T. Voigtmann^, M. E. Cates^ and M. Fuchsi 
^Fachbereich Physik, Universitat Konstanz, D-78457 Konstanz, Germany 
^SUPA, School of Physics, The University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, 



UK 



We consider the nonlinear rheology of dense colloidal suspensions under a time-dependent simple 
shear flow. Starting from the Smoluchowski equation for interacting Brownian particles advected by 
shearing (ignoring fluctuations in fluid velocity) we develop a formalism which enables the calculation 
of time-dependent, far-from-equilibrium averages. Taking shear-stress as an example we derive 
exactly a generalized Green-Kubo relation, and an equation of motion for the transient density 
correlator, involving a three-time memory function. Mode coupling approximations give a closed 
constitutive equation yielding the time-dependent stress for arbitrary shear rate history. We solve 
this equation numerically for the special case of a hard sphere glass subject to step-strain. 

PACS numbers: 82.70.Dd, 64.70.Pf, 83.60.Df, 83.10.Gr 



The constitutive equations of a material relate its stress 
tensor (T{t) at time t to its flow history, described by a 
strain-rate tensor K{t' < t) . Dense suspensions of colloidal 
particles, close to the glass transition, show strongly vis- 
coclastic behaviour and nonlinear shear response' their 
constitutive equations must be highly nontrivial In 
this paper, we report first-principles work aimed at elu- 
cidating these equations theoretically. 

Three alternative approaches to the rheology of glassy 
materials have recently been considered. The phe- 
nomenological soft glassy rheology (SGR) model |2| mim- 
icks the glass transition by a one-particle hopping dy- 
namics in a well chosen distribution of trap-depths 
While the predictions of SGR are broadly consistent with 
experiments on many non-ergodic soft materials it 
does not capture the discontinuous jump in yield stress 
on glass formation observed in experiments on colloidal 
suspensions of hard spheres 0. The same is true of 
spin glass approaches which describe a different phe- 
nomenology akin to 'power-law yielding materials' 0. 
The observed yield-stress is captured however by a first- 
principles approach to colloid rheology, based on Mode 
Coupling Theory (MCT), which has recently been formu- 
lated for systems under steady shear 0|. The MCT has 
had considerable semiquantitative success in accounting 
for the interaction dependence of the static glass tran- 
sition 0, 111 and the time dependence of light scattering 
correlators from microscopic starting points. MCT 
also gives sensible predictions for the viscoelastic spec- 
trum G* (w) as measured in linear rheology 0, ■ 

In this Letter we develop a first-principles description 
of the far-from-equilibrium states of dense colloids un- 
der prescribed time dependent shear flow. For a system 
of interacting particles advected by the imposed strain 
rate k (this neglects velocity fluctuations and hence hy- 
drodynamic interactions |^), we begin by developing a 
formalism which facilitates calculation of general time 
dependent averages. Using this formalism we derive for- 
mally exact generalized Green-Kubo equations, taking 
the shear stress as an example. These expressions can 
be approximated in terms of transient two-time density 
correlators. We find an exact equation of motion for the 



correlator which, for general strain rate, displays a re- 
markable three-time structure to the memory function. 
Finally, we make an MCT-based closure of this equation. 
To show that the theory yields sensible predictions for 
strongly time-dependent fiows, we consider step-strain as 
a specific example. Our approach is valid for all homoge- 
neous and incompressible flows which satisfy K{t)-K{t)=0 
(this enforces homogeneous states with translationally in- 
variant spatial correlations). For clarity of presentation 
we focus on the case of simple shearing with fixed axes 
(velocity along x, gradient along y), and an arbitrary 
time-dependent strain rate 7(i), so that = i{t)5ix^jy. 

Our findings highlight the formal importance of 'in- 
tegration through transients' (ITT) in preparing the 
best ground for judicious application of MCT. Indeed, 
a somewhat simpler MCT-inspired approach to colloid 
rheology was developed in 0,^3, which for fluid states 
in steady shear gives broadly similar results to those of 
However, its recent {ad-hoc) extension to time depen- 
dent flows [l3j gives results quite different from ours; in 
particular, there is no sign of a third time in the memory 
function which we show to be an exact consequence of the 
Smoluchowski equation and which is preserved in our ap- 
proximations. This feature (also missed in another recent 
simplified MCT approach [Ij|) turns out to be crucial in 
capturing the nonlinear response, even for simple time- 
dependent flows. Our flndings should also help guide 
the development of 'fully schematic' single-mode models. 
These models can addre ss p hysics not considered here, 
such as shear thickening [l5| , so their extension to time- 
dependent flows is highly desirable. In future, such mod- 
els may suggest improvements to the MCT factorizations 
used here, iterating towards a fully predictive theory of 
colloid rheology. 

We start with a system of N spherical Brownian par- 
ticles of diameter d interacting via internal forces F^, 
I = 1, • • • , TV, and dispersed in a solvent with a specified 
time dependent velocity profile v(r, t) = K{t)-r. The dis- 
tribution function of particle positions evolves according 
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to the Smoluchowski equation 0, ^3 
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where flit) is the Smoluchowski operator and we have in- 
troduced dimensionless variables for length, energy and 
time, d = ksT = Dq = 1. Translational invariance of 
the sheared system leads to a coupling between a den- 
sity fluctuation with pre-advected wavevector q(t, t') = 
q — q • J^, ds k(s) at time t' and another with wavevector 
q at later time t. Wavevector advection is how strain 
enters our formalism and accounts for the afiine part of 
the particle motion in the imposed flow. 

Integrating the Smoluchowski equation we obtain the 
following formal solution for the distribution function 



^it) = 



dt' Tr{/t(t)<T} e 
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where a is the zero wavevector limit of the potential 
part of the stress tensor and is the equilibrium dis- 
tribution function. The adjoint Smoluchowski operator 
can be obtained by partial integration and is given by 
(t) = Ei {9z + Fi+r,-K'^ it)] ■ di . The time ordered ex- 
ponential e_ imposes that later times appear on the right 
and arises because the operator flit) does not commute 
with itself at different times. Eq. (|2Jl generalizes the 
ITT method developed in and is to be used with the 
understanding that functions to be averaged are placed 
to the right of the operators and then integrated over 
particle coordinates. A general function / of the phase 
space coordinates thus has the time dependent aver- 

age {f)it) = {f)+f_^dt'{TT{Kit)a}cxp_[j!'dsnHm 
with respect to the distribution function |(2Jl, where (•) in- 
dicates averaging with respect to the equilibrium distri- 
bution. By choosing / = axy/V ^3 obtain an exact 
generalized Green-Kubo relation for the time dependent 
shear stress (in volume V) 



ait) 



dt'^it') 
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where the factor [•] can be identified formally as 
G(t,t', [7]), a time-dependent shear modulus. Replacing 
51^ it) with the quiescent-state operator recovers linear re- 
sponse. Eq. Q opens a route to calculate cr(t) for a given 
flow history 7(t). The ITT method based on Eq. ^ also 
yields expressions for correlators, distorted structure fac- 
tors and susceptibilities which will be detailed elsewhere. 

To approximate our formally exact result (jSJ we now 
project onto densities pq ~ X^i ^''^ and density pairs 
(given by the square of the density in real space) [y, 01 ■ 
This physical approximation amounts to assuming that 
these are the only slow variables, sufficient to describe 
the relaxation of the local structure in the glassy regime. 



The resulting shear stress is given by 



ait)^ dt'^it') 
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where 3>k(i,i') is the transient density correlator which 
(in the absence of time translational invariance) is a 
function of two times, and S'f. — dSk/dk with Sk the 
equilibrium static structure factor. The projection onto 
density pairs means that the interparticle forces are 
fully determined from Sk and density fluctuations. This 
MCT-based approximation is well-tested, although if 
the equal-time structure under shear deviates strongly 
enough from Sk to enter an anharmonic regime, im- 
provements to it may be needed 0, In Eq. Q), 
the term ^iit,t') can be viewed as the 'survival prob- 
ability' to time t of a stress contribution created by an 
initial step-strain applied at earlier time t'; the remain- 
ing factor is the stress per unit initial strain. The tran- 
sient density correlator required in Eq. (0J is defined as 
$q(t,f) = (p* eicp_[J^,dsn1is)]pq^t.t'))/iNSg) and is a 
key quantity within our approach. Note that it only con- 
tains information on the strain accumulated between the 
two correlated times t' , t and is independent of the strain 
history for times earlier than t'. 

Equation gives the stress in terms of <I>k(i, t') which 
is itself dependent on flow history. Using Zwanzig-Mori 
type projection operator manipulations and applying the 
theory of Volterra integral equations, we obtain the fol- 
lowing formally exact results: 



d_ 

di 
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^ dt'mq(t,t',to)^$q(t',to)) =0, 



Tqit,to) = - {p*^nsit,to)p'q) 



mq(i, i',to) = 



S, 

1 + S„ , ,2 
— - q ■ Ktt,, ■ q + |q • Ku„ \ 

-(p*^},(^^ to)uit, t', to)Q,n,it, to)p'^) 
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(7) 

In the equation of motion ||SJ) for ^i^it,t'), the initial 
decay rate Tqit, t') describes Taylor dispersion which en- 
hances diffusion in the direction of flow [T^ . For our 
chosen flow geometry, this can be calculated explicitly 
as given in Eq. (O, where = J^^ ds/«(s) is the shear 
strain accumulated between and t. The memory func- 
tion rriqit^t' jtg) in Eq. |S1 describes competition between 
shearing and the cage effect responsible for slow struc- 
tural relaxation; strikingly, this is a function of three 
times, not two. It is useful to interpret mq(<, t' , to) as de- 
scribing the decay of memory between times t' and t, in 
the presence of shear, allowing for the coupling to stress 
in the system that is still relaxing from the strain accu- 
mulated between tg and t' . The time enters the theory 
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in a parametric fashion and is quite distinct in charac- 
ter from the two later times. Eq. |(7J) for rriq involves 

the propagator U{t,t',tf)) = exp_ J^, ds'Qs^"'^{s' ,to), 
where fl^Ht, t') is the single-particle-irreducible opera- 
tor I^, 01 and Qs is an equilibrium projector orthogo- 
nal to density fluctuations. In deriving these formal re- 
sults we introduced Pq = e^'^''^', the density of a single 
tagged particle, whose motion is described by fls{t, tg) = 
-iq • • - iq • Ku^ • (29^ F^) + fjt (t) - q • ■ K^Ug ' 1- 
To close Equations (|S}{ZIl for the transient density cor- 
relator wc now make an MCT-based approximation to 
the average in Eq. 0. Taking care to preserve transla- 
tional invariance of mq{t,t' ,to), we obtain 




mq{t,t',to) 



5,rq(i',to)rq(t,to) 

KZ = q • [ kc^ + pCp + Sq{qCq ~qcq)] 
^qkp = q • [ kCfc + pCp + Ktt„ ■ qSq{ck + Cp 

-^(cgCfe-f CgCp + CfcCp)}], (8) 



FIG. 1: The long time stress a{oo) as a function of strain 7 
for a hard sphere glass just above the MCT glass transition 
{(j> = 0.52, e = {(j> — (j>c)/4'c = 0.008) following a step-strain. 
The dashed line is the linear response result and the solid 
line is a{oo) calculated using the isotropic approximation I2H . 
The inset shows the decaying stress as a function of time 
following the step for 7 = 0.2, 0.3, 0.4 (full lines, bottom to 
top) and 7 = 0.5, 0.6 (broken hnes,top to bottom). 



where p 



k; p = N/V; pcq = 1 - 1/Sg. The 



wavevectors q = q — q • ^tto ^-nd q = q — q • i^w contain 
shear strains accumulated over different temporal ranges. 

Equations H4I5I6() and (jSj) form a closed set of equations 
to predict the shear stress for arbitrary time-dependent 
shear flows of the form k = 'j{t)6ixSjy. Other than 7(t), 
the only required inputs are density p and the static 
structure factor Sk in the unsheared state. The paramet- 
ric nature of to is made explicit in the MCT approxima- 
tion (jH)) where all three times enter the vertex functions 
]/(i.2) Q^-^Yy t' and t enter the correlators. For 7 = 
our equations reduce to those of quiescent MCT [l^ and 
for steady shear to those of Q . 

In developing nonlinear constitutive equations it is 
helpful to study nonlinear step-strain benchmark. 
In an idealized step-strain experiment the shear rate is 
given by 7(i) = ^8{t — Iq)^ which provides a demand- 
ing test of any constitutive equation. For step-strain our 
approximate Green-Kubo relation reduces to 



<7{t) ^Jdk 
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where t > to and k'*' = {kx,ky — ^^kx,kz)- We have 
included an additional superscript on the correlator to 
make explicit the strain dependence. The initial de- 
cay rate becomes independent of time, F^^ = / Sq + 
<lx(lyl{^ + Sq)/[2Sq) + q^"f'^ as the time drops out in 
favour of an explicit 7 dependence. A similar reduction 
occurs for the memory function, leading us to modify the 
notation, mq(t, i', ^q) "i^(t, t'). The memory function 
is thus given in step-strain by 



ml{t,t') 



dk 



C f7 p7 



with 



' 2 ^^ix^ jy 



hatted variables made unhatted, and 
in the vertex expressions of Eq. ||SJ). 



(10) 



The strain is zero in the two correlators since there is 
no strain imposed between t' and t; these arc quies- 
cent MCT correlators. Thus all the 7 dependence in 
the memory function stems from the presence of to in 
Eq. (jSJ. In Eq. (|10|l the forces represented by the ver- 
tices are strain-dependent but relaxation of the struc- 
ture is 7- independent |23|. We thus obtain a simplified 
equation of motion for the correlators needed in 0. At 
statepoints for which quiescent MCT predicts a glass, 
the memory function remains non-ergodic for all values 
of 7. The equation of motion for ^'^{t,to) is a linear 
equation in ^^{t,to) with given non-Markovian memory 
function. Due to the zero duration of the applied strain 
we can make the replacements $^(t,to) — > ^q{t — to), 
m'^{t,t') m'^{t — t') without further approximation. 

In practice the presence of anisotropy still poses dif- 
ficulties for numerical computation; specializing to the 
case of hard spheres we have therefore solved Eq Q for 
the shear stress a{t) within the isotropic approximation 
IT^ l2TI |. The Sk input is taken from the Percus- 
Yevick (FY) theory, which yields a glass transition at 
^nict = T^cPp/Q = 0.516 with our present numerical dis- 
cretization [23. In Figure n we show the long time stress 
a{oo) as a function of strain amplitude 7 for a state point 
just above the MCT glass transition (0 = 0.52, for fluid 
states wc of course find a{oo) ~ 0.) As 7 is increased we 
first obtain the expected linear response behaviour be- 
fore entering a regime of sub-linear increase, indicating 
the onset of plastic flow. (We find that a regime of linear 
response persists even at the critical point (j) = 4>c-) For 
large 7 values our calculations yield unphysical negative 
(t(oo). This may reflect a shortcoming of MCT, where 
for large 7 use of Sk to proxy the colloidal interactions 
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becomes questionable. The inset shows a{t) for various 
values of 7. In contrast to polymer melts following a 
step-strain, (j{t) for the present hard sphere system is 
not strain factorable 0| . 

The rcsuhs (Fig. P) are in qualitative agreement with 
recent step-strain experiments on suspensions of PMMA 
particles above the glass transition |23 | . The experimen- 
tal data show a peak and a region of negative slope in 
cr(c!o) in accord with our findings The isotropic 

approximation is known to underestimate the effects of 
shearing a fuller treatment might shift the peak in 
cr(oo) from 7 ~ 0.4 closer to the experimental results 
[23, which peak at 7 ~ 0.1. Note also that in the exper- 
iments the strain is ramped up over some short finite in- 
terval, during which additional plastic rearrangment may 
occur. Our general expressions H4I5I6I8|) should capture 
this although (|9I10() clearly do not. 

To summarize, for interacting Brownian particles ad- 
vected by a non- fluctuating shear flow, we have general- 
ized the integration through transients formalism of to 
address arbitrary far-from-equilibrium, time-dependent 
shearing. When complimented with mode-coupling ap- 
proximations this provides a route to calculating time de- 
pendent averages in the sheared system. As a demonstra- 
tion of this approach we have presented a Green-Kubo- 
type relation for the shear stress and an equation of mo- 



tion for the transient density correlator, involving a mem- 
ory function with non-trivial time dependence. MCT clo- 
sure of these expressions yielded a first-principles consti- 
tutive model for the shear rheology of dense suspensions 
close to the glass transition. In step-strain, this predicts 
first a linear regime followed by plastic deformation of 
the glass with a maximum in the long-time stress, as 
seen experimentally. 

The formalism excludes extensional flow, and we 
assumed a spatially homogeneous strain rate j{t). 
Nonetheless, the work goes far beyond linear response 
[Tll| : applying the theory to the case of oscillatory shear 
would predict strain dependent storage and loss moduli 
including 'higher harmonic' contributions 0, Effi- 
cient numerical algorithms to tackle the anisotropy and 
loss of time translation invariance are currently under 
development. Meanwhile the formal developments pre- 
sented here form a secure starting point both for more 
complete theories, in which the MCT assumptions might 
be partially relaxed, and for schematic models that sim- 
plify the algebra but add extra physics such as anhar- 
monicity, shear-thickening and jamming fl5l |. 
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